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Abstract 

In this paper we deal with the packings derived by horo- and hy¬ 
perballs (briefly hyp-hor packings) in the n-dimensional hyperbolic 
spaces H” (re = 2,3) which form a new class of the classical packing 
problems. 

We construct in the 2— and 3—dimensional hyperbolic spaces hyp- 
hor packings that are generated by complete Coxeter tilings of degree 
1 i.e. the fundamental domains of these tilings are simple frustum 
orthoschemes and we determine their densest packing configurations 
and their densities. 

We prove that in the hyperbolic plane (re = 2) the density of the 
above hyp-hor packings arbitrarily approximate the universal upper 
bound of the hypercycle or horocycle packing density 4 and in 
the optimal configuration belongs to the [7,3, 6] Coxeter tiling with 
density ps 0.83267. 

Moreover, we study the hyp-hor packings in truncated orthosche¬ 
mes [p, 3,6] (6 < p < 7, p E R) whose density function is attained its 
maximum for a parameter which lies in the interval [6.05,6.06] and the 
densities for parameters lying in this interval are larger that ~ 0.85397. 
That means that these locally optimal hyp-hor configurations provide 
larger densities that the Boroczky-Florian density upper bound (« 
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0.85328) for ball and horoball packings but these hyp-hor packing 
configurations can not be extended to the entirety of hyperbolic space 

m3. 

1 Introduction 

The packing and covering problems with solely horo- or hyperballs (horo- or 
hypespheres) are intensively investigated in earlier works in the n-dimensional 
[n > 2) hyperbolic space H"'. 

1. On horoball packings 

In the n-dimensional hyperbolic space there are 3 kinds of the ’’balls 
(spheres)” the balls (spheres), horoballs (horospheres) and hyperballs 
(hyperspheres). 

The 2-dimensional case of circle and horocycle packings was settled by 
L. Fejes Toth in [B]. 

The greatest possible density in hyperbolic space is ~ 0.85328 which 
is not realized by packing regular balls. However, it is attained by a 

_3 

horoball packing of El where the ideal centers of horoballs lie on the 

_3 

absolute hgure of El . This ideal regular simplex tiling is given with 
Coxeter-Schlafli symbol [3,3,6] see e.g. [1], |1], |3] and |5]. 

In the previous paper [T3] we proved that the above known optimal 
ball packing arrangement in is not unique. We gave several new ex¬ 
amples of horoball packing arrangements based on totally asymptotic 
Coxeter tilings that yield the Boroczky-Florian packing density upper 
bound |1]. Furthermore, by admitting horoballs of different types at 
each vertex of a totally asymptotic simplex and generalizing the sim- 
plicial density function to H” for (n > 2) we hnd the Boroczky type 
density upper bound is no longer valid for the fully asymptotic simplices 
in cases n > 3 [2T], [22]. For example, the density of such optimal, lo¬ 
cally densest packing is ~ 0.77038 which is larger than the analogous 

—4 

Boroczky type density upper bound of ~ 0.73046 for El . However 
these ball packing conhgurations are only locally optimal and cannot 
be extended to the entirety of the hyperbolic spaces H”. 

In the paper [T5| we have continued our investigation of ball packings in 
hyperbolic 4-space using horoball packings, allowing horoballs of differ- 
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ent types. We have shown seven connterexamples (which are realized 
by allowing one-, two-, or three horoball types) to a conjectnre of L. 
Fejes-Toth abont the densest ball packings in hyperbolic 4-space. The 
maximal density is ~ 0.71645 

In [31] we proved that the optimal horoball density related to the hy¬ 
perbolic 24 cell in is ~ 0.71645 as well. 

2. On hyperball packings 

In [25] and [26] we have stndied the regnlar prism tilings and the cor¬ 
responding optimal hyperball packings in H” (n = 3,4) and in the 
paper [^ we have extended the in former papers developed method 
to 5-dimensional hyperbolic space and constrnct to each investigated 
Coxeter tiling a regnlar prism tiling, have stndied the corresponding op¬ 
timal hyperball packings by congruent hyperballs, moreover, we have 
determined their metric data and their densities. 

In hyperbolic plane the universal upper bound of the hypercycle 
packing density is - proved by I. Vermes in [33] and recently, (to the 
author’s best knowledge) the candidates for the densest hyperball (hy¬ 
persphere) packings in the 3,4 and 5-dimensional hyperbolic space H" 
are derived by the regular prism tilings which are studied in papers 

[ 25 ], [ 25 ] and [27] . 

In the universal lower bound of the hypercycle covering density is 
determined by I. Vermes in [34] . 

In the paper [28] we have studied the n-dimensional (n > 3) hyperbolic 
regular prism honeycombs and the corresponding coverings by con¬ 
gruent hyperballs and we have determined their least dense covering 
densities. Moreover, we have formulated a conjecture for the candidate 
of the least dense hyperball covering by congruent hyperballs in the 3- 
and 5-dimensional hyperbolic space. 

In [30] we studied the problem of hyperball (hypersphere) packings in 
the 3-dimensional hyperbolic space. We described to each saturated hy¬ 
perball packing a procedure to get a decomposition of the 3-dimensional 
hyperbolic space into truncated tetrahedra. Therefore, in order to 
get a density upper bound to hyperball packings it is sufficient to de¬ 
termine the density upper bound of hyperball packings in truncated 
simplices. We considered the hyperball packings in truncated sim- 
plices and proved that if the truncated tetrahedron is regular, then 
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the density of the densest packing is ~ 0.86338 which is larger than the 
Boroczky-Florian density npper bound, however these hyperball pack¬ 
ing conhgurations are only locally optimal and cannot be extended to 
the entirety of the hyperbolic spaces 

In this paper we deal with the packings with horo- and hyperballs (briefly 
hyp-hor packings) in the n-dimensional hyperbolic spaces H"' (n = 2, 3) which 
form a new class of the classical packing problems. 

We construct in the 2— and 3—dimensional hyperbolic spaces hyp-hor 
packings that are generated by complete Coxeter tilings of degree 1 i.e. the 
fundamental domains of these tilings are simple frustum orthoschemes with 
a principal vertex lying on the absolute quadric Q and the other principal 
vertex is outer point. We determine their densest packing conhgurations and 
their densities. These Coxeter tilings exist in the 2—, 3— and 5—dimensional 
hyperbolic spaces (see [10]) and have given by their Coxeter-Schlahi graph 


in Fig. 1. 
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the hyperbolic plane 
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Figure 1: Coxeter-Schlahi graph of Coxeter tilings of degree 1. 


above hyp-hor packings arbitrarily approximate the universal upper bound 
of the hypercycle or horocycle packing density ^ and in the optimal 
conhguration belongs to the [7, 3, 6] Coxeter tiling with density ~ 0.83267. 

Moreover, we consider the hyp-hor packings in truncated orthoschemes 
[p, 3, 6] (6 < p < 7, p G R). Its density function is attained its maximum 
for a parameter which lies in the interval [6.05, 6.06] and the densities for 
parameters lying in this interval are larger that ~ 0.85397. That means that 
these locally optimal hyp-hor conhgurations provide larger densities that the 
Boroczky-Florian density upper bound (~ 0.85328) for ball and horoball 
packings but these hyp-hor packing conhgurations can not be extended to 
the entirety of hyperbolic space 
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2 The projective model and 
the complete orthoschemes 

For H"' we use the projective model in the Lorentz space of signature 
( 1 , n), i.e. denotes the real vector space equipped with the bilinear 
form of signature (1, n) ( x, y) = —x^y^ + + ■ • • + x'^y'^ where the non¬ 
zero vectors x = (a;°, , x"') G and y = {y^, y^,..., y'^) G 

are determined up to real factors, for representing points of P”(M). Then HI"’ 
can be interpreted as the interior of the quadric Q = {[x] G P”|( x, x) = 
0 } =: 9EI” in the real projective space Vn+i). 

The points of the boundary 5IHI"' in P” are called points at inhnity of H"', 
the points lying outside dW^ are said to be outer points of relative to Q. 
Let P([x]) G P”, a point [y] G P” is said to be conjugate to [x] relative to 
Q if { X, y) = 0 holds. The set of all points which are conjugate to P([x]) 
form a projective (polar) hyperplane pol{P) := {[y] G P”|( x, y) = 0}. Thus 
the quadric Q induces a bijection (linear polarity —)■ Vn+i)) from the 

points of P"' onto its hyperplanes. 

The point V[x] and the hyperplane a[a] are called incident if xa = 0 
(xGV-+i\{0}, aG V„+i\{0}). 

Definition 2.1 An orthoscheme S in H"' (2 < n G N) is a simplex bounded 
by n -I- 1 hyperplanes ..., such that (see [21 [11]) H'-PW, for j ^ 
i — + 1. 

The orthoschemes of degree d in H” are bounded by n -|- d -|-1 hyperplanes 
..., such that id* T for j 7 ^ i — 1 , i -f 1 , where, for 

d = 2, indices are taken modulo n -|- 3. For a usual orthoscheme we denote 
the (n -|- l)-hyperface opposite to the vertex Aj by id* (0 < i < n). An 
orthoscheme S has n dihedral angles which are not right angles. Let a*-^ 
denote the dihedral angle of S between the faces id* and HP Then we 
have a*-^ = |, if 0 < i < j — 1 < n. The n remaining dihedral angles 
Q,*,*-i-i, (^0 < i < n — 1 ) are called the essential angles of S. Geometrically, 
complete orthoschemes of degree d can be described as follows: 

1. For d = 0, they coincide with the class of classical orthoschemes intro¬ 
duced by Schlafli (see Dehnitions 2.1). The initial and hnal vertices, 
Aq and An of the orthogonal edge-path AjAi+i, i = 0 ,... ,n — 1 , are 
called principal vertices of the orthoscheme. 
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2. A complete orthoscheme of degree d = 1 can be interpreted as an 
orthoscheme with one onter principal vertex, say An, which is trnn- 
cated by its polar plane pol{An) (see Fig. 1 and 3). In this case the 
orthoscheme is called simply truncated with outer vertex An- 

3. A complete orthoscheme of degree d = 2 can be interpreted as an 
orthoscheme with two outer principal vertices, Aq, An, which is trun¬ 
cated by its polar hyperplanes pol{Ao) and pol{An). In this case the 
orthoscheme is called doubly truncated. We distinguish two different 
types of orthoschemes but I will not enter into the details (see [?]). 

A n-dimensional tiling V (or solid tessellation, honeycomb) is an inhnite set of 
congruent polyhedra (polytopes) that £t together to £11 all space (H" (n ^ 2 )) 
exactly once, so that every face of each polyhedron (polytope) belongs to 
another polyhedron as well. At present the cells are congruent orthoschemes. 
A tiling with orthoschemes exists if and only if each dihedral angle of a tile 
is submultiple of 27r (in the hyperbolic plane the zero angle is also possible). 

Another approach to describing tilings involves the analysis of their sym¬ 
metry groups. If V is such a simplex tiling, then any motion taking one 
cell into another maps the entire tiling onto itself. The symmetry group of 
this tiling is denoted by SymV. Therefore the simplex is a fundamental do¬ 
main of the group SymV generated by rehections in its (n — l)-dimensional 
hyperfaces. 

The scheme of an orthoscheme S' is a weighted graph (characterizing 
S C up to congruence) in which the nodes, numbered by 0,1,...,n 
correspond to the bounding hyperplanes of S. Two nodes are joined by an 
edge if the corresponding hyperplanes are not orthogonal. For the schemes of 

k, k2 k„, k„ 

O -O-O. o-o-o 

complete Coxeter orthoschemes S C H"' we adopt the usual conventions and 
sometimes even use them in the Coxeter case: If two nodes are related by 
the weight cos ^ then they are joined by a (p — 2)-fold line for p = 3, 4 and 
by a single line marked p for p > 5. In the hyperbolic case if two bounding 
hyperplanes of S are parallel, then the corresponding nodes are joined by a 
line marked cxo. If they are divergent then their nodes are joined by a dotted 
line. 

The ordered set [ki,..., fcn-i, ^n] is said to be the Coxeter-Schlahi symbol 
of the simplex tiling V generated by S. To every scheme there is a corre- 
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spending symmetric matrix (c*-^) of size (n + 1 ) x (n + 1 ) where c** = 1 and, 
for z ^ j G {0,1, 2,..., n}, eqnals — cos ^ with all angles between the 
facets i,j of S. 

For example, (c*-^) below is the so called Coxeter-Schlafli matrix of the 
orthoscheme S in 3-dimensional hyperbolic space with parameters (nodes) 
ki=p,k2 = q,h = r 



/I — cos ^0 ^ \ 

— cos - 1 — cos - 0 

P Q 

0 - cos I 1 - cos f 

V 0 0 — cos ^ 1 / 


(2.4) 


3 Basic notions and formulas 

3.1 Coxeter tilings generated by simply frustum or¬ 
thoschemes 

In general the complete Coxeter orthoschemes were classihed by Im Hof in 
|9] by generalizing the method of Coxeter and Bohm, who showed that they 
exist only for dimensions < 9. From this classiheation it follows, that the 
complete orthoschemes of degree d = 1 exist up to 5 dimensions. 

In this paper we consider the orthoschemes of degree 1 where the initial 
vertex Aq lies on the absolute quadric Q. These orthoschemes and the cor¬ 
responding Coxeter tilings exist in the 2-, 3— and 5—dimensional hyperbolic 
spaces and are characterized by their Coxeter-Schlafli symbols and graphs 
(see Fig. 1). 

In n-dimensional hyperbolic space H"' (n > 2) it can be seen that if 
S = A0A1A2 ■■■ An P0P1P2 ■.. Pn is a complete orthoscheme with degree 
d = 1 (a simply frustum orthoscheme) where An is a outer vertex of H"' 
then the points Pq, Pi, P 2 ,..., Pn-i he on the polar hyperplane vr of An (see 
Fig. 2 in H^). We consider the images of S under reflections on its side 
facets. The union of these n-dimensional orthoschames (having the common 
TT hyperplane) forms an inhnite polyhedron denoted by Q. Q and its images 
under reflections on its ,,cover facets” £11 hyperbolic space H"' without overlap 
and generate n-dimensional tilings T. 

The constant k = \J~^ 'is the natural length unit in H”. K will be the 
constant negative sectional curvature. In the following we assume that k = 1. 
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a. b. 

Figure 2 : a. A 3-dimensional complete ortlioscheme of degree d = 1 (simple 
frustum orthoscheme) with outer vertex A 3 . This orthoscheme is truncated 
by its polar plane tt = pol{A^). b. Two congruent adjacent simple frustum 
orthoschemes. 

3.2 Volumes of the n-dimensional 
Coxeter orthoschemes 

1 . 2 -dimensional hyperbolic space 

In the hyperbolic plane a simple frustum orthoscheme is a Lambert 
quadrilateral with exactly three right angles and its fourth angle is 
acute ^ (g > 3) (see Fig. 1). In our case the Lambert quadrilateral has 
a vertex at the inhnity i.e. the angle at this vertex is 0. Its area can 
be determined by the well-known defect formula of hyperbolic triangles 
(see [7j): 

Voh{S) = (3.1) 

2. 3-dimensional hyperbolic space 

Our polyhedron A 0 A 1 A 2 P 0 -P 1 -P 2 is a simple frustum orthoscheme with 
outer vertex A 3 (see Fig. 1) whose volume can be calculated by the 
following theorem of R. Kellerhals |11] : 

Theorem 3.1 The volume of a three-dimensional hyperbolic complete 
orthoscheme (except Lambert cube cases) S is expressed with the es¬ 
sential angles Ooi) ^ 12 , 0 ^ 23 , (0 < Op < |) (Fig. 1) in the following 
form: 
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^ 0 / 3 ( 5 ) — -{£(q;oi + 6) — C{aQi — 9) + >C(— + 0:12 — 9)+ 

+ ~ *^12 ~ 9) + £( 0:23 + 9) — C{a2z — 9) + — 9)}, (3.2) 

where 9 G [0, |) is defined by the following formula: 


tan(6*) 




cos^ ai2 


sin^ aoi sin^ 0:23 


cos Ooi cos 023 


and where C{x) := — flog \2smt\dt denotes the Lobachevsky funetion. 

0 

For our prism tilings Tpqr we have: Ooi = 012 = ^23 = 7 • 

3.3 On hyperballs 

The equidistant surface (or hypersphere) is a quadratic surface that lies at 
a constant distance from a plane in both halfspaces. The inhnite body of 
the hypersphere is called a hyperball. The n-dimensional half-hypersphere 
{n = 2, 3) with distance h to a hyperplane tt is denoted by Pfif. The volume 
of a bounded hyperball piece Pfifi^An-i) bounded by an (n — l)-polytope 
An-i C TT, Pfif and by hyperplanes orthogonal to tt derived from the facets 
of An-i can be determined by the formulas (3.3) and (3.4) that follow from 
the suitable extension of the classical method of J. Bolyai: 

VohififiiAi)) = VolfiAi) sinh (h), (3.3) 

Vokin’^iA 2 )) = ^V"o/ 2 (. 42 ) [sinh (2h) + 2h] , (3.4) 

where the volume of the hyperbolic (n — l)-polytope An-i lying in the plane 
TT is V oln-l{An-l). 

3.4 On horoballs 

A horosphere in H” [n > 2) is a hyperbolic n-sphere with inhnite radius 
centered at an ideal point on dIP. Equivalently, a horosphere is an (n — 1)- 
surface orthogonal to the set of parallel straight lines passing through a point 
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of the absolute quadratic surface. A horoball is a horosphere together with 
its interior. 

We consider the usual Beltrami-Cayley-Klein ball model of H" centered 
at 0 ( 1 , 0 , 0 , ..., 0) with a given vector basis e* (z = 0 , 1 , 2 ,..., n) and set 
an arbitrary point at inhnity to lie at Tq = (1, 0,..., 0,1). The equation 
of a horosphere with center Tq = (1, 0,..., 1) passing through point S = 
( 1 , 0 ,..., s) is derived from the equation of the the absolute sphere —-f- 
+ • • • + = 0 , and the plane — x" = 0 tangent to the 

absolute sphere at Tq. The general equation of the horosphere is in projective 
coordinates (s ^ ± 1 ): 


(s — 1) I —x°x° -I- y (x*)^ 1 — (1 -I- s)(x^ — x"'Y = 0, 


E 

i=l 


(3.5) 


and in cartesian coordinates setting hi = ^ it becomes 


2(Er=>?), 


+ 


{i-sy 


= 1 . 


(3.6) 


In n-dimensional hyperbolic space any two horoballs are congruent in the 
classical sense. However, it is often useful to distinguish between certain 
horoballs of a packing. We use the notion of horoball type with respect to 
the packing as introduced in [ 22 ] . 

In order to compute volumes of horoball pieces, we use Janos Bolyai’s 
classical formulas from the mid 19-th century: 


1. The hyperbolic length L{x) of a horospheric arc that belongs to a chord 
segment of length x is 

L(x) = 2sinh (3.7) 


2. The intrinsic geometry of a horosphere is Euclidean, so the {n — 1)- 
dimensional volume ^ of a polyhedron A on the surface of the horo¬ 
sphere can be calculated as in The volume of the horoball piece 

'H(A) determined by A and the aggregate of axes drawn from A to the 
center of the horoball is 


-A 


n 


Vol{n{A)) 


(3.8) 
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4 Hyp-hor packings in hyperbolic plane 

We consider the previously described 2-dimensional Coxeter tilings given 
by the Coxeter symbol [cxd] (see Fig. 1), which are denoted by Ta- The 
fundamental domain of 7^ is a Lambert quadrilateral AqAiPiPq (see Fig. 3) 
that is denoted by iFa- It is derived by the truncation of the orthoscheme 
710^1^2 by the polar line vr of vertex A 2 where the initial principal vertex of 
the orthoschemes Aq is lying on the absolute quadric Q and its other principal 
vertex A 2 is an outer point of the model. 

Its images under reflections on its sides £11 hyperbolic plane without 
overlap. The tilings Ta contain a free parameter 0 < a < 1, a G R. The 
polar straight line of A 2 is tt and tt fl AqA 2 = Pq, tJ" C ^1^42 = Pi. 

We consider the usual Beltrami-Cayley-Klein ball model of centered at 
0(1, 0, 0) with a given vector basis {i = 0,1, 2) and set the above Lambert 
quadrilateral AqAiPiPq in this coordinate system with coordinates 

Ao(l,0,l); Lli(l,0,0); Pi(l,a,0); Po(l, a, 1 - a^); (0 < a < 1). 

The polar line ^ 2 ( 1 , — b, 0)^ of the outer vertex ^ 2 ( 1 , b, 0) is tt which contains 
the points Pq and Pi (see Fig. 3). 

We construct hyp-hor packings to Ta tilings therefore the hyper- and 
horocycles have to satisfy the following requirements: 

1. The centre of the horocycle can only be the vertex Aq and the corre¬ 
sponding horocycle SjaiVi) has not common points with inner of seg¬ 
ments AiPi and PqPi. These horocycle types depend on parameter 
a of the considered tiling Ta and passing through the point T(l,0,|/i) 
(0 < a, Hi < 1) (see Fig. 4). 

2. We can choose the base straight line of the hypercycle 'Ha{y 2 ) between 
the lines PqPi and AiPi, the role of these lines is symmetrical regarding 
the packings. We consider the AiPi line as base line of hypercycles 
to construct hyp-hor packings. Furthermore, 'Ha(l/ 2 ) has not common 
points with inner of segments AqPq. These hypercycle types depend 
on the parameter parameter a of the considered tiling Ta and passing 
through the points T(l, 0 , 1 / 2 ) (0 < a, y 2 < 1) (see Fig. 4). 

3. card{int{S^a{yi)) A int{PLa{y 2 ))] = 0. 
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If the hyper-and horocycles hold the above requirements then we obtain hyp- 
hor packings Ta in the hyperbolic plane derived by the structure of the con¬ 
sidered Coxeter simplex tilings. 


A, 



Figure 3: The fundamental domain of 7^ is a Lambert quadrilateral 
AqAiPqPi, at present a = 0.5. 


Definition 4.1 The density of the above hyp-hor packings PaiyiiUi) is 

_ yol{Pa n (ila(2/l) U Ua{yi)) 

b\Pa\y\-,y2)) — -r———-. 

Vol{Pa) 

It is well known that a packing is locally optimal (i.e. its density is locally 
maximal), then it is locally stable i.e. each ball is fixed by the other ones 
so that no ball of packing can be moved alone without overlapping another 
ball of the given ball packing or by other requirements of the corresponding 
tiling. Therefore, we can assume that the horocycle S)a{yi) and the hypercy¬ 
cle 'Ha{y 2 ) touch each other at the point T{l,0,y) where (0 < y < 1). The 
possible values of y = yi = y 2 may depend on parameter a (see Fig. 4). 

4.1 Main types of hyp-hor packings 

We distinguish two main types of hyp-hor packings: 
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a. b. 

Figure 4: a. The hyp-hor packing of type 1, at present a = 0.7. b. The 
hyp-hor packing of type 2, at present a = 0.4. 


1. The hypercycle PLaiy) contains the point Pq and the horocycle P}a{.y) 
touches it at T{l,0,y). These conhguration can be realized for all 
possible parameters 0 < a < 1. These packings are denoted by Vaiu) 
(see Fig 4.a). 

2. The horocycle Sja{y) passes through the point Pq and the hypercycle 
Piaiy) touches it at T{l,0,y). These conhgurations exist if 0 < a < ^. 
If a = then y = 0 i.e. iOa(O) touches the line AiPi at Ai, (the height 
of the hypersphere is 0). If a > then these conhgurations do not 
satisfy the requirements of the hyp-hor packings. These packings are 
denoted by Va{y) (see Fig. 4.b). 


4.1.1 The densities of packings Vliy) 


In this case the coordinates of touching point T{l,0,y) can be easily ex¬ 
pressed as the function of parameter a\ y = \/l — a?. We obtain by the 
formulas (3.1), (3.3), (3.6), (3.7), (3.8) and by Dehnition 4.1 that the density 
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of the packings {y) of type 1 can be calculated by the following formula: 

= 6{V\a)) = 

/ /I , /I 2-2Vl-a^-a^ + 2a^\\ n -, 

( 4 sinh ( - arccosh ( -- - -j j + 2 Vl — a? j /tt 

^ (4.1) 

where 0 < a < 1. 


Lemma 4.2 Analysing the above density formula we obtain that 


lim 

d — ^0 


5{V\a)) 


—, lim 

TT 


5{V\a)) 


2 

TT 


and ^ < 6(V^{a)) < | for parameters 0 < a < 1 (see Fig. 5.a). 


Corollary 4.3 In the hyperbolic plane the universal upper bound density 
of ball packings can be arbitrarily accurate approximate with the densities 
6(V^{a)) of hyp-hor packings of type 1. 


4.1.2 The densities of packings V^iy) 


Similarly to the previous section the coordinates of touching point T{l,0,y) 
can be expressed as the function of parameter a: y = 1 — 2a^. We obtain 
by the formulas (3.1), (3.3), (3.6), (3.7), (3.8) and by Dehnition 4.1 that the 
density of the packings V^iy) of type 2 can be calculated by the following 
formula: 


4 sinh - arccosh — 


5{Vl{x)) = 5{V\a)) = 
1 -3 + 20 ^' 


2 1 -a2 


-1 + 2a^ 

Vl — 


/tt 


(4.2) 


where 0 < a < 

Lemma 4.4 Analysing the above density formula we obtain that 


lim 

CL —^0 L 


5{V\a)) 


lim 

1 


TT L 


5{r\a)) 


2^2 

TT 


and < 5{V^{a)) < | for parameters 0 < a < ^ (see Fig. 5.b). 

Corollary 4.5 In the hyperbolic plane the universal upper bound density 
of ball packings can be arbitrarily accurate approximate with the densities 
5{V‘^{a)) of hyp-hor packings of type 2. 
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Figure 5: a. The density function of hyp-hor packings of type 1. b. 
The density function of hyp-hor packings Pf of type 2. 


4.2 The general cases 

1 . First we consider the hyp-hor packings Va{y) where the conhgurations 
are ’’between the two main cases”: i.e. the inequalities 0 < a < ^ and 

1 — 2a^ < y < \/l — hold. 

2. We get the second case if the inequalities < a < 1 and 0 < ?/ < 
\/l — a'^ hold. 

In both cases the densities of packings Va{y) are denoted by S(Va{y)) which 
can be determined by the formulas (3.1), (3.3), (3.6), (3.7), (3.8) and by 
Dehnition 4.1: 


1 -l + 2|/-2a2 + 2aV-|/^ 

2 a2 (1 - 1/2) 



5{Va{y)) = (^4 sinh ( - arccosh 



(4.3) 

Analyzing the above density function we obtain that the maximal densities 


can be attained at the ’’main cases” described in subsections 4.1.1 and 4.1.2. 
Therefore we get the following 
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Theorem 4.6 In the hyperbolic plane the densest packing configurations 
with horo- and hyperballs generated by simple frustum orthoschemes with 
Schldfly symbol [cxd] provide the Vfiy) {i = 1,2) packings (described in sub¬ 
section f.1.1 and f.1.2) if their parameter a —)■ 0. Their densities SiVfia)) 
are arbitrarily accurate approximate the universal upper bound density | of 
ball packings ofUf. 

Remark 4.7 If y = 0 and < a < 1 then we obtain ball packings Pa(0) 
which contain purely horocycles (see Fig. 6.a). Their densities can be com¬ 
puted also by the formula (4-3) and its graph is illustrated in Fig. 7. The 
maximal density is 0.90032 belonging to a = ^. 




a. b. 

Figure 6: a. The horocycle packing Va{0) where a = ■^. b. The density 
function of horocycle packings Pa(0). 


5 Hyp-hor packings in hyperbolic space 

In the 3-dimensional hyperbolic space there are 3 inhnite series of the simple 
frustum Coxeter orthoschemes with vertex at the infinity listed in Fig. 1 and 
characterized in Sections 2-3. The considered tilings with Schlafli symbol 
[p, q, r], (g, r) = (3, 6), (4,4), (6, 3) are denoted by where p is an integer 
parameter and p > 7 if (g, r) = (3, 6), p > 5 if (g, r) = (4,4), p > 4 if (g, r) = 
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(6,3). The fundamental domain of is a simple frustum orthoscheme 

AQA 1 A 2 P 0 P 1 P 2 (see Fig. 2) that is denoted by It is derived by the 

truncation of the orthoscheme AQA 1 A 2 A 2 , by the polar hyperplane tt of vertex 
A^ where the initial principal vertex of the orthoschemes Aq is lying on the 
absolute quadric Q and its other principal vertex A 3 is an outer point of the 
model. 

Its images under reflections on its faces £11 hyperbolic space without 
overlap. The polar plane of A 3 is n and n fl AiAs = Pi, {i = 0,1, 2) (see 
Fig. 2). 

We consider the usual Beltrami-Cayley-Klein ball model of centered at 
0(1, 0, 0, 0) with a given vector basis e* {i = 0, 1, 2, 3) and set the above sim¬ 
ple frustum orthoscheme in this usual coordinate system (see Fig. 7.a- 

b). 

Po(l, 0,0,0); Pi(l,0,|/,0); P2{l,x,y,0) where 



and Ao(l, 0,0,1); Ai{l,H,y, zi)-, A 2 {l,x,y, Z 2 ), where 


the Ath coordinates of the points Ai and A 2 can be derived by the following 
procedure described in general for n-dimensional hyperbolic space H"': 


1. The points Pk[Pk] and Ak[aik] {k = 1, 2) are proper points of hyperbolic 
n-space and P^ lies on the polar hyperplane po/(A„)[a”] of the outer 
point An thus 


Pfc ~ c • a„ -I- Hfc G a” c ■ Ana^ + = 0 c = 


Pfc ~ ~ a./j(a„(2 ) an(afc(2 ) — a.khnn ^nhkni 

3.r) 


where is the inverse of the Coxeter-Schlafli matrix (e.g. see (2.4) 
in H^) of the considered orthoscheme. 


2 . 


The hyperbolic distance PkAk can be calculated by the following for¬ 
mula: 

cosh PkAk = cosh h = (Pfc; k) _ _ 

V (Pfc) Pfc) (a-fc, a/j) 



(5.3) 



18 


Jeno Szirmai 


3. The coordinates Zk {k = 1,2) can be derived by the following equations 
(see Fig. 7.a-b): 


coshPi^i = 




13 


— 1 + 


COshP2^2 = 


^22 ^33 ~ ^23 


^11 ^33 a /(—1 + a; 2)(—1 + + zj) ’ 

—1 + 


^22 ^33 a /(—1 + 1 + + 1 /^ + 2 :|) 

(5.4) 

For example for tiling the above coordinates are the following (see 

Fig. 7.a-b): 

1 x/S-i/108 cos^ (y ) — 81 1 \/3a/— 3 + 4 cos^ (4) 

X =--PS 0.27580, y = --Ri 0.47770, 

18 cos(f) ’ ^ 2 cos(f) 

TT 


= 12 \/2 J-566 cos2 j + 318 + 157 cos j cos^ (^yj ps 0.92394, 


^2 = 12 ( — 1 + 2 cos"^ 


— cos 


18 - 31 cos2 


+ 8 cos 


0.69574. 


We construct hyp-hor packings to 7^*''^’^^ tilings therefore the hyper- and 
horocycles have to satisfy the following requirements: 


1. The centre of the horoball can only be the vertex Aq and the corre¬ 

sponding horoball (horosphere) ^p^’^\ 7 i(p)) has not common points 
with inner of side faces P 0 P 1 P 2 and P 1 P 2 A 1 A 2 of the simple frustum 
orthoschem = P 0 P 1 P 2 A 0 A 1 A 2 (see Fig. 7.a-b). These horoball 

types depend on the metric data of the considered tiling and it is 
passing through the point Ti(l, 0, 0, 7 i(p)) where (see (5.3) and (5.4)) 
0 < 7i(p) < 1- 

2. The base plane of the hyperball (hypersphere) 'k 6 p’''\ 72 {p)) is P 1 P 2 P 3 

plane. Furthermore, the above hyperball has not common points with 
inner of side faces AqAiA 2 . These hyperball types depend on the metric 
data of the considered tiling and it is passing through the point 

T 2 ( 1 , 0, 0, 72(p)) where 0 < 71 (p) < tanh(arcosh(P 2 ^ 2 )) < 1 (see (5.3) 
and (5.4)) (in some cases the points Ti and T 2 can coincide i.e. the 
corresponding horoball and hyperball touch each other). 
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a. b. 

Figure 7: The simple frustum orthoscheme with Schlafli symbol [7,3,6] in 
Beltrami-Cayley-Klein model. 


3. card{int{S)^p’'^\'yi{p))) n int{'Hp’''\'y 2 {p)))} = 0. 

If the hyper-and horocycles hold the above requirements then we obtain hyp- 
hor packings Tp^’^^(7i(p), 72 (p)) in the hyperbolic plane derived by the struc¬ 
ture of the considered Coxeter simplex tilings. 

Definition 5.1 The density of the above hyp-hor packings (7i(p)) 72 (p)) 

is 




Voi{4‘^’^^ n {4p‘''\ii{p))un4\y2{p))) 




5.1 Hyp-hor packings to and tilings 

5.1.1 On tilings 

It is well known that a packing is locally optimal (i.e. its density is locally 
maximal), then it is locally stable i.e. each ball is hxed by the other ones 
so that no ball of packing can be moved alone without overlapping another 
ball of the given ball packing or by other requirements of the corresponding 
tiling. Therefore, hrst we consider the largest possible horo- and hyperballs 
to the considered tiling. 
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The largest possible horoball centered at Aq is passing through the vertex 
Ai and the largest possible hyperball contains the vertex A 2 . We get by easy 
calculations, that these ’’maximal large balls” have common inner points 
for any permissible parameters p. Thus, the optimal arrangement can be 
achieved if the horo- and hyperballs touch each other i.e. holds the z(p) = 
7i(p) = l 2 {p) equation and its density of the hyp-hor packing Vp'^’'^\z{p)) 
depends on parameter z{p) where 0 < z{p) < tanh(arcosh(P2^2)) < 1 (see 
((5.3) and (5.4)). The volumes of VA S^^p’'^^){z{p)) and Vol{J^p^’'^^ fl 
'h!'p'^'^){z{p)) can be calculated by the formulas (3.2), (3.4), (3.6), (3.7), (3.8), 
(5.3), (5.4) for any parameters R 3 p > 4 but only to parameters Z 3 p > 5 
belong tilings and so hyp-hor packings in The largest possible horoball is 
denoted by S)!'p’'^\zmax{p)) passes through the point Ti = T2(l, 0, 0, 2:max(p)) 
which is the common point with the hypersphere 'h}p''^\zmax{p))■ We blow 
up this hypersphere keeping the horoball S)i^p''^\zmax{p)) tangent to it upto 
this hypersphere touches the faces A1A2A3 at vertex A 2 . During this process 
we can compute the densities of considered packings as the function of 0 < 
C(p) < arcosh(P2^2)) — sirtanh.{zmax{p)) (see (5.4)) by the following equation 




n (ai,*'‘>(C(p)) u (C(p))) 


artanh(2;(p)) = aitanh^Zmax {p)) + C{p)- 


For example, if p = 5 then 0 < C(5) < tanh(arcosh(P2A2)) — Zmax{5) ~ 
0.33419 and the graph of (5(P5^’^^(C(5))) is described in Fig. 8.a. Analyzing 
the above density function we get that the maximal density is achieved at the 
endpoint of the above interval with density 0.81296 (see Table 1). Similarly 
to the above p = 5 parameter we can compute the optimal densities for all 
possible other parameters p. The results for some parameters are summarized 
in Table 1. 

The densest hyp-hor configuration among the considered ’’realizable” pack¬ 
ings belongs to packing with density zz 0.812956. 


Table 1, q = A, r = A 

P 




5 

0.34084197 

0.27709010 

0.81295769 

6 

0.38165233 

0.30003810 

0.78615556 

7 

0.40369221 

0.30777518 

0.76240058 
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a. b. 

Figure 8: The graphs of density functions and (5(Pg^’^^(C(6))). 

5.1.2 On tilings 

The determination of the optimal hyp-hor packing configurations of packings 
7^(6,3) (^2 3 p > 4) is similar to the above tilings therefore here we only 
summarize the results in Table 2. 

The densest hyp-hor configuration among the considered packings belongs 
to packing with density ^ 0.81209. 


Table 2, g = 6, r = 3 

p 



6(1^) 

4 

0.31716925 

0.25756985 

0.81208961 

5 

0.35991902 

0.27187731 

0.75538469 

6 

0.38060310 

0.27009741 

0.70965634 


5.2 Hyp-hor packings to tilings 

The investigations of these tilings are a little different from the above tilings. 

First we consider the largest possible horo- and hyperballs to the consid¬ 
ered tiling. 

The largest possible horoball centered at Aq is passing through the vertex 
Ai and the largest possible hyperball contains the vertex A 2 . Contrary to 
the above tilings here we get by easy calculations, that these ’’maximal 
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large balls” have not common points for any permissible parameters p. 

The volumes of V hi Vol{J'p^ fl )(72(p)) 

can be calculated by the formulas (3.2), (3.4), (3.6), (3.7), (3.8), (5.3), (5.4) 
for any parameters R 3 p > 6. 

Remark 5.2 VTe note here, that only the parameters Z 3 p > 7 provide 
tilings and hyp-hor ball paeking configurations in the hyperbolic space 
The other parameters provide locally optimal density as well but these hyp- 
hor packing configurations can not be extended to the entirety of hyperbolic 
space . 

We can compute the densities for all possible parameters p of hyp-hor pack¬ 
ings. The results for some parameters are summarized in Table 1. 

The densest hyp-hor configuration among the considered ’’realizable” pack¬ 
ings belongs to packing with density ~ 0.83267. 


Table 3, g = 3, r = 6 

p 



Jil^) 

7 

0.31781165 

0.26463185 

0.83266882 

8 

0.34695830 

0.27901923 

0.80418664 

9 

0.36482363 

0.28351212 

0.77712105 


Finally, we obtain by careful computations and investigations from the 
above method and results the following 

Theorem 5.3 The packing configuration (see Section 5.2) provides 

the maximal density ~ 0.83267 of hyp-hor packings ({q,r) = (4,4), 

(6, 3), (3, 6) andp is suitable integer parameter (see Fig. 1)) which are derived 
by the Coxeter tilings generated by complete orthoschemes of degree 1 (simple 
frustum orthoschemes). 

5.2.1 On non-extendable hyp-hor packings Vp (6 < p < 7, p G R) 

The computation method described in the former sections is suitable to de¬ 
termine the densities of hyp-hor packings for parameters (6 < p < 7, p G R) 
as well. To any parameter belongs a simple frustum orthoscheme and there¬ 
fore we can determine similarly to the above cases the corresponding density 
of its optimal hyp-hor packing. But these packings can not be extended 
to the 3-dimensional space. Analyzing these non-extendable packings for 
parameters (6 < p < 7, p G R) we obtain the following 
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Figure 9: The optimal hyp-hor packing configuration to Coxeter tiling 
[7,3,6], 


Theorem 5.4 The function (6 < p < 7, p G R) is attained its 

maximum for the parameter Popt which lies in the interval [6.05, 6.06] and 
the densities for parameters lying in this interval are larger that ~ 0.85397. 
That means that these locally optimal hyp-hor configurations provide larger 
densities that the Bdrdczky-Florian density upper hound (~ 0.85328) for ball 
and horoball packings (l^). 

Remark 5.5 We note here, that the b-dimensional analogous periodic hyp- 
hor packing will be investigated in a forthcoming paper. 

The question of finding the densest hyp-hor packing without any symme¬ 
try assumption in the n-dimensional hyperbolic space is open. Similarly to 
it, the discussion of the densest horoball and hyperball packings in the n- 
dimensional hyperbolic space n > 3 with horoballs of different types and 
congruent hyperballs has not been settled yet (see [H], [15], [21], [22]). 

Moreover, optimal sphere packings in other homogeneous Thurston ge¬ 
ometries represent another huge class of open mathematical problems. For 
these non-Euclidean geometries only very few results are known (e.g. [T9] , 
ra, m, PD- By the above these we can say that the revisited Kepler 
problem keep several interesting open questions. 





24 


Jeno Szirmai 


References 

[1] Bezdek, K. Sphere Packings Revisited, Eur. J. Combin., 27 jQ 
(2006), 864-883. 

[2] Bohm, J - Hertel,E. Polyedergeometrie in n-dimensionalen 
Rdumen konstanter Kriimmung, Birkhauser, Basel (1981). 

[3] Boroczky, K. Packing of spheres in spaces of constant curvature, 
Acta Math. Acad. Sci. Hungar., 32 (1978), 243-261. 

[4] Boroczky, K. - Florian, A. Uber die dichteste Kugelpackung 
im hyperbolischen Raum, Acta Math. Acad. Sci. Hungar., 15 
(1964), 237-245. 

[5] Fejes Toth, G. - Kuperberg, G. - Kuperberg, W. Highly Satu¬ 
rated Packings and Reduced Goverings, Monatsh. Math., 125/2 
(1998), 127-145. 

[6] Fejes Toth, L. Regular Figures, Macmillian (New York), 1964. 

[7] G. Horvath, A. Formulas on hyperbolic volume, Aeguat. Math., 
83/1 (2012), 97-116. 

[8] Johnson, N.W., Kellerhals, R., Ratcliffe, J.G., Tschants, S.T. 
The Size of a Hyperbolic Goxeter Simplex, Transform. Groups, 
4/4 (1999), 329-353. 

[9] Im Hof, H.-G. A class of hyperbolic Goxeter groups. Expo. Math., 
(1985) 3 , 179-186. 

[10] Im Hof, H.-G. Napier cycles and hyperbolic Goxeter groups. 
Bull. Soc. Math. Belgigue, (1990) 42 , 523-545. 

[11] Kellerhals, R. On the volume of hyperbolic polyhedra. Math. 
Ann., (1989) 245 , 541-569. 

[12] Kellerhals, R. Ball packings in spaces of constant curvature and 
the simplicial density function, J. Reine Angew. Math., (1998) 
494,189-203. 



Packings with horo- and hyperballs ... 


25 


[13] Jacquemet, M. The inradius of a hyperbolic truncated n- 
simplex, Discrete Comput. Geom., 51/4 (2014), 997-1016 DOI; 
DOI 10.1007/s00454-014-9600-y. 

[14] Kozma, R.T., Szirmai, J. Optimally dense packings for fully 
asymptotic Coxeter tilings by horoballs of different types, 
Monatsh. Math., 168/1 (2012), 27-47. 

[15] Kozma, R.T., Szirmai, J. New Lower Bound for the Optimal 
Ball Packing Density of Hyperbolic 4-space, Discrete Comput. 
Geom., 53/1 (2015), 182-198, DOI; 10.1007/s00454-014-9634-l. 

[16] Molnar, E. The Projective Interpretation of the eight 3- 
dimensional homogeneous geometries, Beitr. Algebra Geom.,, 
38/2 (1997), 261-288. 

[17] Rogers, C.A. Packing and Covering, Cambridge Tracts in Math¬ 
ematics and Mathematical Physics 54, Cambridge University 
Press, (1964). 

[18] Szirmai, J. The optimal ball and horoball packings to the 
Coxeter honeycombs in the hyperbolic d-space, Beitr. Algebra 
Geom., 48/1 (2007), 35-47. 

[19] Szirmai, J. The densest geodesic ball packing by a type of Nil 
lattices, Beitr. Algebra Geom., 48/2 (2007), 383-397. 

[20] Szirmai, J. The densest translation ball packing by fundamental 
lattices in Sol space, Beitr. Algebra Geom., 51/2 (2010), 353- 
373. 

[21] Szirmai, J. Horoball packings to the totally asymptotic regular 
simplex in the hyperbolic n-space, Aeguat. Math., 85 (2013), 
471-482, DOI; 10.1007/s00010-012-0158-6. 

[22] Szirmai, J. Horoball packings and their densities by generalized 
simplicial density function in the hyperbolic space, Acta Math. 
Hung., 136/1-2 (2012), 39-55, DOI; 10.1007/sl0474-012-0205- 



26 


Jeno Szirmai 


[23] Szirmai, J. Regular prism tilings in SL(2,R) space, Aequat. 
Math., (2014) 88/1-2, 67-79, DOI; 10.1007/s00010-013-0221-y. 

[24] Szirmai, J. Simply transitive geodesic ball packings to x R 
space groups generated by glide reflections, Ann. Mat. Pur. 
Appl, (2014) 193/4, 1201-1211, DOI; 10.1007/sl0231-013- 
0324-z. 

[25] Szirmai, J. The p-gonal prism tilings and their optimal hyper¬ 
sphere packings in the hyperbolic 3-space, Acta Math. Hungar. 
(2006) 111 (1-2), 65-76. 

[26] Szirmai, J. The regular prism tilings and their optimal hyperball 
packings in the hyperbolic n-space, Publ. Math. Debrecen (2006) 
69 (1-2), 195-207. 

[27] Szirmai, J. The optimal hyperball packings related to the 
smallest compact arithmetic 5-orbifolds, Submitted Manuscript 
(2013). 

[28] Szirmai, J. The least dense hyperball covering to the regular 
prism tilings in the hyperbolic n-space, Ann. Mat. Pur. Appl. 
(2014), DOI; 10.1007/sl0231-014-0460-0. 

[29] Szirmai, J. A candidate for the densest packing with eqnal balls 
in Thnrston geometries, Beitr. Algebra Geom. (2014) 55/2, 441- 
452, DOI; 10.1007/sl3366-013-0158-2. 

[30] Szirmai, J. Hyperball packings in hyperbolic 3-space, Submitted 
Manuscript (2014). 

[31] Szirmai, J. Horoball packings related to hyperbolic 24 cell. Sub¬ 
mitted Manuscript (2015). 

[32] Vermes, I. Uber die Parkettiernngsmoglichkeit des dreidimen- 
sionalen hyperbolischen Raumes durch kongrnente Polyeder, 
Studia Sci. Math. Hungar. (1972) 7, 267-278. 

[33] Vermes, I. Ausfiillnngen der hyperbolischen Ebene durch kon- 
gruente Hyperzykelbereiche, Period. Math. Hungar. (1979) 
10/4, 217-229. 



Packings with horo- and hyperballs ... 


27 


[34] Vermes, I. Uber regulare Uberdeckungen der Bolyai- 
Lobatschewskischen Ebene durch kongruente Hyperzykel- 
bereiche, Period. Math. Hungar. (1981) 25/3, 249-261. 

Budapest University of Technology and Economics Institute of Mathematics, 

Department of Geometry, 

H-1521 Budapest, Hungary. 

E-mail: szirmai@math.bme.hu 
http://www.math.bme.hu/ '^szirmai 



